Correlation amplitude for S = | XXZ spin chain in the critical region: 
Numerical renormalization-group study of an open chain 
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The density-matrix renormalization-group technique is used to calculate the spin correlation 
functions {SjS%} and {S^S^) of the one- dimensional S = i XXZ model in the gapless regime. 
The numerical results for open chains of 200 spins are analyzed by comparing them with correlation 
functions calculated from a low-energy field theory. This gives precise estimates of the amplitudes of 
the correlation functions in the thermodynamic limit. The exact amplitude recently conjectured by 
Lukyanov and Zamolodchikov and the logarithmic correction in the Heisenberg model are confirmed 
numerically. 
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The quantum spin chains have been a subject of very 
active research over the many years. Among others, the 
one-dimensional S = h XXZ model, 
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is a simplest nontrivial model and has served as a test 
ground of various theoretical techniques. These include 
the exact calculation via the Bethe ansatz,tl numerical 
computations such as the exact diagonaUzation studies,EI 
and the conformal field theory (CFT).H Although each 
approach alone is often not powerful enough to provide 
sufficient information we need, employing these tech- 
niques together has proven to be very successful to study, 
e.g., finite-size energy spectra and spin correlation func- 
tions. It is well known that the equal-time correlation 
functions of the XXZ chain show the power-law decay 
in the critical regime (—1 < A < 1) at zero tempera- 
ture. In the thermodynamic limit [L — > oo) they have 
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where only the leading oscillating and non-oscillating 
terms are written. The exact expression of the parameter 
r\ appearing in the exponents, 



r\ = 1 cos 1 A, 
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was obtained by comparing the result of the effective- 
field-theory description-, (abelian bosonization) with the 
Bethe ansate solutionis On the other hand, little had 
been knowntl on the correlatiou.amplitudes, until recently 
Lukyanov and Zamolodchikovl (LZ) conjectured the ex- 
act form of A x : 
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where £ = r)/(l — f]). This result was then used 
in Refs. |^,^ to obtain the correlation amplitude for 
the antiferromagnetic Heisenberg model (A = 1) with 
the logarithmic correction predicted earlier from the 
renormalization-group argument .0il3 The aim of this pa- 
per is to numerically determine the correlation ampli- 
tudes (A x , Ax, and A z ) for broad ranges of the anisotropy 
A in the critical regime. To this end we apply the density- 
matrix renormalization-group (DMRG) methodEJ-o to 
compute the spin correlation functions for spin chains 
of L = 200, and fit the numerical data to functional 



forms expected from the effective field theory (abelian 
bosonization). Indeed, the idea of fitting the correlation 
functions of finite periodic systems to the CFT form has 
been applied successfully to the Heisepbipc case (A = 1) 
to examine the logarithmic correction.liaiLZlln the present 
study we use open spin chains instead of periodic ones, 
not only because of the better performance of the DMRG 
method for open chains, but because the boundary ef- 
fects provide useful information. We numerically verify 
the LZ formula (^) and give numerical data for the yet 
analytically unknown amplitudes A x and A z . 
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We consider a finite XXZ chain of L spins with open 
boundaries, Eq. ([l]). We assume L to be an even in- 
teger so that the ground state is singlet. According to 
the standard abelian bosonization method,Ert3 the low- 
energy dynamics of the XXZ model is described by the 
Gaussian model 
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where <p(x) and </>(x) are bosonic fields. We identify the 
site index j with the continuous variable x, and impose 
boundary conditions 0(0) = 4>(L + 1) = const, to respect 
the fact that there is no spin at j = 0, L + l. The fields 
have the mode expansions 
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where q n 
and [a rn , 
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This leads to the commutation 
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relation [0(x),0(y)] = -(i/2)[l + sgn(a: - j/)]. The 
ground state |0) is the singlet vacuum state of the bosons: 
a n|0) = Q|0) = 0. The spin operators in the original 
Hamiltonian ([!]) can be expressed in terms of the phase 
fields as 
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where a. b, and c are real constants. It follows that 



S? = c(-l) j cos[27ri?0(a;)] - ib sin[27r R(j>{x)\ sin • 
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Note that the second term in Eq. j7c[ ) is hermitian be- 
cause [<fi(x) , (fi(x)] = —i/2. These formulas are slightly 
modified from those in Refs. |l8| |2(] for reasons which will 
become clear below. It is then straightforward to eval- 
uate the correlation functions for the ground state |0), 
yielding 
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and (gtj) we used the following regu- 
- cos(q„:r)]/n = ln[/i(x)]. Note that 
reduce to Eqs. (2a) and ( ^b| ) with 
2 /i a „A J - in tlie limit 

<I. It is also im- 
) and (^tj) have exactly 



In deriving Eqs 
larization: 
Eqs. (fa|) and" 

A z = a 2 : /2, A x = c 2 /2, and .1 
L — > 00 with |j - || <L and 
portant to observe that Eqs. 

the same j- and fc-dependence at the Heisenberg point 
(rj = 1) in accordance with the SU(2) spin symmetry. In 
fact, under the given boundary conditions, it is essential 
to use the boson representation (|7a|) and j7q) to have this 
SU(2) property. Equation (8b) has another nice feature 
that, when r\ = 1/2 and a — I/71", it coincides with the 

exact correlation function (SjS%) of the XY model. We correlation amplitudes, and the results are summarized 
use Eqs. (|8a|) and (|B) to analyze the numerical data. 



We calculated the spin-spin correlation functions 
(SfS%) (a = x, z) for the ground state of L = 200 
spin chains using the DMRG method with the improved 
algorithmOiiJ We employed the finite system method, 
and the maximum number of kept states m is 80. From 
the difference between the data at m = 80 and m — 60, 
we estimate the numerical error due to the truncation of 
the Hilbert space to be of order 10~ 5 for (5J Sf) and 10~ 6 
for (5JiS|). Figures [j] and || show some of our numeri- 
cal data of the correlations between and S%, where 
j = (L - r + l)/2 and k = (L + r + l)/2 for odd r and 
j = (L- r) 12 and k = (L + r)/2 for even r. The DMRG 
results are shown by open symbols, whose size is, how- 
ever, made larger than the numerical error for illustra- 
tion. Taking a, &, and c as free parameters and using the 
exact exponent Eq. (||), we fit the datacil using Eq. ( |8a|) 
and (p9)- The small dots connected by lines in Figs. Ill 
and gare the results of the fitting. From the relation 
72, A x = 6 2 /4, and A z = a 2 /2, we obtain the 
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in Table I. The excellent agreement between the LZ's 
conjecture and the value obtained from the fitting 
clearly demonstrates that Eq. (^) is indeed exact, and 
that there is little finite-size correction left in a, b, and 
c because Eqs. ( |8a| ) and (3b) correctly account for most 
of the finite-size effect. We also note that the numerical 
estimate of A z at A = agrees with the exact value, 
1/27T 2 . For A < —0.7, it is difficult to estimate A x and 
A z because the exponent I/77 for these terms becomes 
so large that they give only a negligible contribution to 
the correlation functions. As seen in Figs, [l] and |^, for 
— 1 < A < 0.6 Eqs. ( |aj) and Jsb| ) can fit the numerical 
data extremely well including the oscillations of period 
4, which are characteristic of open chains. 



10 



10 



A 

to _ 
-10 

to 
V 



10 



10" 




200 



10" 



A 
W 

to" 



10" 



10 




200 



FIG. 1. (-iy- k {S*S£) versus r = \j - k\ for A = -0.5, 0, 
and 0.5 (from above). The open symbols are the DMRG data 
and the small dots connected by lines are the fits, Eq. (|8a 



FIG. 2. \{SjSk) I versus r = \j - k\ for A = 0.5 (above) and 
—0.5. The open symbols are the DMRG data and the small 
dots connected by lines are the fits, Eq. (|8q). 
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where g is determined from 
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with 7 = 0.5772 . . .. Equation (|10|) is shown as a 
solid line in Fig. |[ which goes on top of the numerical 
data for r ^ 50. Thus, our numerical result is consis- 
tent with the recent analytic calculation,Q'El (S"S%) = 

[log(2 % /2^e 7+1 r)/87r 3 ] 1 / 2 /r for r > 1. The discrepancy 
at larger r in Fig. || would be due to the finite-size ef- 
fect or the boundary effect. The numerical estimates of 
the amplitudes for A = 0.8 and 0.9 also suffer from the 
correction coming from the irrelevant operator and have 
larger error bars than for A < 0.6. 



As A becomes closer to 1, however, the quality of the 
fitting becomes poor. This is due to the irrelevant op- 
erator cos(2</>/i?) ignored in Eq. (||). Roughly speaking, 
this operator gives rise to corrections of order r -4 ' 1 /'' -1 ' 
to the amplitudes, andj-at A = 1 this turns into the 
logarithmic correctionE'El Since this type of correction 
is not taken into account in Eqs. ( |8a| ) and (|t]), they 
cannot describe the correlation functions completely for 
the parameter regime where the effect of the leading ir- 
relevant operator becomes important. To elucidate this, 
we plot in Fig. || the effective amplitude A z (r) given by 
A z (r) = a 2 (r)/2, where a(r) is obtained by solving Eq. 
( |Sh| ) for each r = \j — k\. We see that A z (r) is con- 
stant for A = 0.5, which is another indication that the 
fitting works well. At A = 1 , on the other hand, A z (r) is 
roughly an increasing function of r. This result is com- 
pared with the recent analytical result valid in the limit 
L — > 00 ,@ 



0.1 












A=l 


a 




— field theory 


ocP 






A=0.5 


oca 








1 







50 100 150 200 

r 

FIG. 3. Effective correlation amplitude A z (r) as a function 
of r for A = 1 and 0.5. The solid line represents the y/log r 
correction [Eq. (^)] predicted from a field-theory calculation. 
A z (r) is almost constant (~ 0.1) at A = 0.5. 
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In summary, we have demonstrated that the spin cor- 
relation functions of open XXZ chains in the critical 
regime can be fitted very well with the correlation func- 
tions calculated from the effective low-energy theory (free 
bosons), if the effect of the irrelevant operator is not 
severe. From the fitting, we obtained the correlation 
amplitudes for broad ranges of the parameter A with 
high precision and numerically verified the Lukyanov and 
Zamolodchikov's conjecture. We hope that our numeri- 
cal estimates for A z and A x will be confirmed in turn by 
analytic calculations in the future. The numerical com- 
putation was carried out at the Yukawa Institute Com- 
puter Facility, Kyoto University. A.F. is supported by a 
Monbusho Grant for overseas research. 
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TABLE I. Correlation amplitudes estimated from the fitting of the numerical data. A x is the conjectured amplitude, 
Eq. (M). The figures in parentheses indicate the error bar on the last quoted digits. 



A 


six 




Ax 


A z 


-0.9 


0.15567 


0.15560(3) 






-0.8 


0.15904 


0.1589(1) 






-0.7 


0.15968 


0.1595(1) 




0.008(1) 


-0.6 


0.15912 


0.1589(2) 


0.0048(14) 


0.0133(1) 


-0.5 


0.15786 


0.1576(2) 


0.0076(9) 


0.0184(4) 


-0.4 


0.15617 


0.1558(2) 


0.0099(7) 


0.0235(2) 


-0.3 


0.15417 


0.1538(3) 


0.0122(4) 


0.02921(3) 


-0.2 


0.15196 


0.1515(3) 


0.0144(2) 


0.03556(3) 


-0.1 


0.14958 


0.1491(3) 


0.0164(2) 


0.0425(2) 


0.0 


0.14709 


0.1466(4) 


0.0182(2) 


0.0501(5) 


0.1 


0.14451 


0.1440(4) 


0.01995(7) 


0.0588(3) 


0.2 


0.14187 


0.1414(3) 


0.02154(5) 


0.0683(6) 


0.3 


0.13921 


0.1388(3) 


0.02296(4) 


0.0791(8) 


0.4 


0.13656 


0.1362(3) 


0.02420(4) 


0.0918(9) 


0.5 


0.13400 


0.1339(1) 


0.02525(6) 


0.1063(9) 


0.6 


0.13164 


0.1320(1) 


0.0261(2) 


0.1236(5) 


0.7 


0.12973 


0.1308(5) 


0.0267(3) 


0.145(1) 


0.8 


0.12896 


0.132(1) 


0.0271(7) 


0.171(5) 


0.9 


0.13214 


0.138(3) 


0.027(2) 


0.20(1) 
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